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ABSTRACT 
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We made optical observations of the Gegenschein using a liquid-nitrogen-cooled 
wide-field camera, Wide-field Imager of Zodiacal light with ARray Detector (WIZARD), 
between March 2003 and November 2006. We found a narrow brightness enhancement 
superimposed on the smooth gradient of the Gegenschein at the exact position of the 
antisolar point. Whereas the Gegenschein morphology changed according to the orbital 
motion of the Earth, the maximum brightness coincided with the antisolar direction 
throughout the year. We compared the observed morphology of the Gegenschein with 
those of models in which the spatial density of the interplanetary dust cloud was consid- 
ered and found that the volume scattering phase function had a narrow backscattering 
enhancement. The morphology was reproducible with a spatial distribution model for 
infrared zodiacal emission. It is likely that the zero-phase peak (the so-called opposi- 
tion effect) was caused by coherent backscattering and/or shadow- hiding effects on the 
rough surfaces of individual dust particles. These results suggest that big particles are 
responsible for both zodiacal light and zodiacal emission. Finally, we derived the geo- 
metric albedo of the smooth component of interplanetary dust, assuming big particles, 
and obtained a geometric albedo of 0.06±0.01. The derived albedo is in accordance 
with collected dark micrometeorites and observed cometary dust particles. We con- 
cluded that chondritic particles are dominant near Earth space, supporting the recent 
theoretical study by dynamical simulation. 

Subject headings: interplanetary medium — instrumentation: miscellaneous — comets: 
general — minor planets, asteroids 



1. Introduction 



The Gegenschein is a faint glow in the sky in the antisolar direction. Several hypotheses had 
been suggested to explain the Gegenschein: sunlight reflected by particles concentrated near the 
libration point in the Sun-Earth-particles system, emission from the Earth's gaseous tail, and 
sunlight scattered by the Earth's dust tai l. After the detectio i i of the Gege nschein at 1.86 AU 
from the Sun by the Pioneer 10 spacecraft (jWeinberg et al.lll973l : Illannerlll974l ). it has been widely 
accepted that the Gegenschein is a portion of zodiacal light enhanced by the backscattering of 
interplanetary dust particles. Intensive observations of the Gegenschein date back nearly half a 
century. It was observed by integrating the brig htness iri the f ield of view using photomultipliers 



attached to telescopes o r telephotographic lenses (jTanabe 



1975 



1965 



Roosenlll970al : iDumont &: Sanchez 



Kwon et al.l l2004l ) . These authors observed the sky background within a certain interval 



on the sky (l°-5°) or scanned the sky with large-field-of-view instruments. A volume scattering 
phase function of the interplanetary dust parti c les was deduced based on low-spatia l -resolution data 



(IDumont &: Sanchejll975l : iLeinert et al.lll976l : iHong 
19861 ). 



1985 



Lumme &: Bowell 



1985 



Lamv &: Perrin 
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Modern obser vation techniques enable us to obtain snapshot images of the Gegenschein. 
James et al.l (jl997l ) were the first to used a cooled CCD camera attached to a fish-eye lens for 
the observation of zodiacal light and the Gegenschein, obtaining sufficient resolution (4.8' pixel"^) 
images. Subsequently, a research group at Kobe University obtained images using a Peltier-cooled 
CCD camera (with a quantum efficiency of 10% at 0.44 fim) wi th an equatorial mou nting and 



1999 



30). 



Note, 



succeeded in detecting asteroidal dust bands on the Gegenschein (jlshiguro et al. 
however, that they applied 1° x 1° boxcar smoothing to reduce the pixel-to-pixel noise level, which 
could blur image resolution in many cases. The recent space observation with the Solar Mass Ejec- 
tion Imager (SMEI) prov ided nearly all-sky pho tometric m aps for over 5 years wit h a moderate 
spatial resolution of 30' (jBuffington et al.l |2009|) . Note that iBuffington et al.l ((20091) stacked data 
over the mission period because the S/N ratio was inadequate to produce the Gegenschein map. 
They did not consider the seasonal variation in the Gegenschein morphology. 

The first question we would like to address in thi s paper i s the ex act position of the maximum 
brightness in the Gegenschein. Dumont (jl965l ) and Roosen ( 1970a b) foun d that the maximum 
bright ness was close to the antisolar point. However, Tanabel ~( 1965 ) and Wolstencroft Sz Rose 
(j 19671 ) reported that the maximum brightness of the Gegenschein was displaced frorn the antisolar 
point and existed close to the invariable plane of the Solar System. iMukai et al.l ()2003l ) argued 
that the maximum brightness was largely influenced by the existence of fine-scale structures in th e 



zodiacal cloud, i.e., asteroidal dust bands (jSvkes fc Greenbere]ll986l : ISvkeslll990l : lReach et al.lll997l ) 



based on their CCD observation. After subtraction of the fine structures, they esti mated a deviation 
from the antisolar point of —0.4° from an image taken on 2 November 1997. IBuffington et al 



(j2009l ) found that the maximum existed at the antisolar point, although the study was based on 
S MEI measurenients, whos e spatial resolution ( 0.5°) is equivalent to the displacement discussed 



m 



Mukai et al.l (120031 ) and llshiguro et al.l (jl998l ). High-resolution observations (<0.5°) are thus 



essential to settle the matter of the maximum brightness position. 

The second question c oncerns t he G egenschein morphology. A systematic morphological 
change was pointed out bv iTanabd (119651^ and later con firmed by high-resolution observations 
(jjames et al.l 119971 : llshiguro et al.l Il998l : iKwon et al.l l2004l ) . This change can be caused by the 



discrepancy between the orbital plane of the Earth (i.e., the ecliptic plane) and the plane of sym- 
metry of the zodiacal cloud. The spa tial distribution of the zodiacal cloud has been investigated 
around the solar elongation e ~ 90° (iKelsall et al.lll998l ). However, little is known about the spa- 
tial distribution around the Gegenschein region, which might retain information about the spatial 
distribution of interplanetary dust particles that exist further than those observed around e ~90°. 

In this paper, we present high-resolution images of the Gegenschein taken by a wide-field 
CCD camera for the diffuse night sky. Wide-field Imager of Zodiacal light with ARray Detector 
(WIZARD). Bot h the sensitivity and the spatial resolution are better than those of the instruments 
used previouslv (James et al.lll997l : llshiguro et al.lll998l : iMukai et al.ll2003l : IBuffington et al.ll2009l : 
Kwon et al.ll2004l ) . We examined the morphology and the photometric center of the Gegenschein for 
4 years using WIZARD. We succeeded in detecting the narrow brightness enhancement at the exact 
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position of the antisolar direction. The maximum brightness of the Gegenschein consistently existed 
at the antisolar point throughout our observation period. We examined the volume scattering phase 
function around the opposition and considered the reason for the narrow brightness enhancement. 
In Section 2, we describe the observations and the data reduction. In Section 3, we show the 
morphological maps of the Gegenschein. In addition, we discuss the annual variation in the peak 
position. In Section 4, we test several phenomenological three-dimensional models to interpret the 
observational results. Finally, we discuss the geometric albedo of interplanetary dust particles in 
the last section. 



2. Observations 
2.1. The Instrument and its Performance 

All data in thi s pape r were taken wi t h WI ZAR D. A detailed descr iption of the instrument is 



given m lUsui et al.l (j2002l ) . Ilshiguro et al.l (j2003l ). and lUeno et al.l (|2007l ). In this section, we briefly 
describe the instrument and report on its performance based on the data reduction results. 

WIZARD was developed to cover a wide-extended zodiacal light and the Gegenschein. It 
consists of a camera head and a wide- field lens (Figured]). The system was mounted on an equatorial 
mounting. We employed a back-illuminated full-frame CCD chip, EEV CCD 42-80 (27.6 x 55.3 
mm; 2048x4096 pixels), placed in a compact N2-cooled cryogenic dewar (Kadel Engineering Corp.). 
During observations, the CCD temperature was maintained at around 180 K. The optical system 
was designed by the Genesia Corporation. It has a large image circle (60 mm in diameter) to cover 
the whole imaging area of the CCD. The focal ratio and the focal length are 2.8 and 32.5 mm, 
respectively. With the camera head, it gives a field of view of 46° x 92° with a pixel resolution 
of 1.4'. Figure [2] shows the transmittance of the optical filter used for zodiacal light observations. 
It covers a wavelength range of 440-520 nm, avoiding prominent airglow emissions and artificial 
emission light. In particular, it was designed not to transmit O I airglow at 557.7 nm and artificial 
mercury vapor light at 435.8 nm. In addition, sodium beacons at 589.2 nm as laser guide stars are 
nowadays problematic for astronomical observations. The transmittances at these wavelengths are 
negligible, that is, 0.35% at 435.8 nm (Hg), 0.19% at 557.7 nm O I, and 0.0008% at 589.2 nm (Na), 
respectively. The CCD clock pattern and the operation commands were prepared on a Linux host 
computer and sent to the CCD chip through a COGITO-3 system. COGITO-3 is a programmable 
control system for imaging devices specifically designed for scientific purposes and was developed 
by Ueno and Wada at ISAS/JAXA. It contains three units: a sequencer board, an A/D board, and 
a driver board. The obtained analog signal was converted to the digital one on the COGITO-3 
A/D board and transferred to the host computer via a PCI-VME bus interface adapter. 

The analog signal in each pixel was sampled 16 times to reduce the readout noise. In addition, 
during CCD readout, we turned off the pulse motor used for star-tracking to reduce noise associated 
with the motor drive unit. Consequently, the readout noise was reduced to the level of 20-25 e~. 
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To monitor the zero level of the CCD, we recorded the overscan data by reading 50 dummy pixels 
before and after each row was read. For confirmation, we examined the stability of the zero intensity 
level using dark frames obtained at the intervals of the light frame. We found that the standard 
deviation of the zero level in each frame was <10 e~. Using the overscan region, we determined 
the zero level to an accuracy of 4 e~ . The typical brightness of the Gegenschein plus airglow was 
about 1.1 e~ s~^. With an exposure time of 600 s, the uncertainty of the zero intensity level (i.e., 
4 e~) was small enough (0.6% of the Gegenschein plus airglow brightness, i.e., 660 e~) to derive 
the brightness of the Gegenschein plus airglow. In addition, background noise becomes dominated 
by photon noise from the natural night sky background rather than the readout noise. For this 
reason, we fixed the exposure time of 600 s for our Gegenschein observations. 

An image taken by a wide-field camera generally shows strong vignetting. The detected in- 
tensity around the optical center is thus brighter than that at the edge of the frame. A flat-field 
correction is essential to obtain accurate intensity maps of the zodiacal light and the Gegenschein. 
It is, however, difficult to obtain flat data using WIZARD because of its huge field of view. We 
investigated flat-fielding for different spatial scales, that is, by considering optical vignetting as a 
low-frequency component and pixel-to-pixel variation in the sensitivity as a high-frequency com- 
ponent. We examined the optical vignetting by observing bright stars at different positions on the 
CCD chip. We confirmed that the weather conditions were stable and the airmasses of these stars 
were constant during the exposures, suggesting that the brightnesses of these light sources were 
stable enough to examine the vignetting. Figure [3] shows the vignetting function across the CCD 
chip. From the measurement, we found that the detected intensity was almost constant within 20° 
from the optical center and reduced by 15% at 38° and 30% at 47°. For reference, we compare 



the yignetting function wit h those of previous CCD observations of the Gegenschein (jjames et al 



19971 : llshiguro et al.lll999al ). in which commercial products for single- lens reflex cameras were used. 
It is clear that the optics for WIZARD is well designed to suppress the vignetting throughout a 
large portion of the CCD chip. Small-scale (pixel-to-pixel) fluctuations of the CCD sensitivity were 
examined by taking images with a Labsphere integrating sphere. We found that the pixel-to-pixel 
responsivity variation was negligibly small (<0.3%). 

The performance of WIZARD, which are critical for estimating the error in the measurements, 
are summarized in Table 1. 



2.2. Observations and Data Analysis 

In normal operation, we began observation soon after the end of evening twilight. We took 
a couple of images successively without slewing the camera to its next position. Since the motor 
controller for star-tracking was turned off for about 3 min after exposures for the CCD readout, the 
position of celestial objects moved about one degree on the CCD chip. This observation strategy 
enables us to distinguish small-scale structures of the Gegenschein from instrumental artifacts 
associated with uneven sensitivity of the instrument. 
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We made a sequence of observations of the zodiacal light and the Gegenschein with WIZARD 
between March 2003 and November 2006 on Mauna Kea. Our observations were supported by 
the Subaru telescope (National Astronomical Observatory of Japan) and the InfraRed Telescope 
Facility (ITRF; NASA and Institute for Astronomy, University of Hawaii). We observed zodiacal 
light and the Gegenschein for 45 nights under photometric conditions. We selected Gegenschein 
data taken at high Galactic latitude near the zenith. In total, 14 nights of data in different seasons 
of the year were analyzed for the present study. We summarize the observation conditions of these 
data in Table 2. Gegenschein data were not available in May-July and December-January because 
the antisolar point was close to the Galactic plane (Figure H]). 

The night sky brightness observed from the ground consists of those of zodiacal light includ- 
ing the Gegenschei n, integrated starlig ht of unresolved stars, diffuse Galactic light, and pseudo- 
continuum airglow (jLeinert et al.lll998l ). These light sources are attenuated by scattering of atmo- 
spheric molecules and aerosols. At the same t ime, the scattered component in Earth's atmosphere 
results in a diffuse light source. iDumontI ()l965l ) introduced the concept of ef fective optic a l dept h for 
diffuse light sourc es, which is different from that for point sources. Later, iHong et al.l (jl998l ) and 
Hong et al.l (j2002l ) developed the idea by solving radiative transfer in an anisotropically scattering 
spherical atmosphere of the Earth. It is difficult reduce zodiacal light data in high-airmass (i.e., 
large zenith distance) regions owing to the considerable contributions of the atmospheric diffuse 
light and the strong zenith extinction, although some intensive studies were co nducted to obtain 



the absolute b rightness of the zodiacal light from ground-based observatories (jKwon et al 



2004 



Dumontlll965l ). The data reduction of the Gegenschein is, however, relatively easy. The Gegen- 
schein appears near the zenith around midnight on Mauna Kea, where the contaminations from 
airglow and zenith extinction are smallest in the sky. The typical brightness of atmospheric diffuse 
light components was 20 Sqiq and had a weak slope of 10%-20% of the Gegenschein brightness 
over a 40° field of view in Gegenschein exposures, where 1 S01O i s equivalent to the f lux of a solar- 
type star of magnitude 10 distributed over one degree squared ([Leinert et al.lll998l ). Inadequate 
reduction for these atmospheric components does not tend to lead to a misunderstanding of the 
morphology of the Gegenschein. In addition to these diffuse light sources mentioned above, there 
remain a number of point sources. Stars brighter than a visual magnitude of 12-14 (depending 
on the Galactic latitude) were resolved by our obser vations. We subtra cted point sources brighter 
than magnitude 12 using the USNO-B2.0 catalog (IMonet et al.ll2003l ). The detected pixel data 
were replaced by the median values of circumjacent pixels. 

For the most part, the data were analyzed as follows: First, the raw data were reduced using the 
overscan region as a bias level plus flat-field data. The fiux calibration factor and the optical depth of 
the point sources are determined by aperture photometry of field stars at various airmasses. Since we 
applied a nonstandard filter, we constructed a photometric standard star catalog for WIZARD data. 
An empirical formula is applied to convert from the Johnson-Cousins system into the WIZARD 
system, that is, MwzD = V + a{B — V), where a is the color conversion factor and B and V are 
B-band and V-band magnitudes, respectively. To derive the coefficient a, we selected 56 stars in 
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BSC Ver. 5 in an airmass range of 1.00 to 1.07. Since the typical extinction coefficient on Mauna 
Kea in the WIZARD band was ~0.15, atmospheric extinction did not change the stellar magnitude 
by ^0.01 magnitude in this sky area, which is adequate to derive the coefficient a. The result is 
shown in Figure [5l In the figure, the bottom axis denotes the instrumental magnitude measured by 
aperture photometry, and the perpendicular axis denotes the magnitude in the WIZARD system 
{Mwzd)- The best-fit parameter was found to be a = 0.605 it 0.005. 

Unresolved starlight brightness plus diffuse Galactic light brightness were subtracted by us- 
ing Pioneer 10 and 11 data, which was obtaine d beyond the inne r Solar System (i.e., outside 
the zodiacal cloud) and the USNO-B2.0 catalog (|Monet et al.ll2003l ). Although the central wave- 
length of WIZARD (460 nm) is slightly longer than that of the Pioneer blue band (440 nm), 
we used the Pioneer data without color correction. The contribution from unresolved starlight 
(stars fainter than magnitude 12) was estimated to be 20 Sqio (or 10% of the Gegenschein bright- 
ness) at a Galactic latitude of 25° and 10 Sqiq (or 5% of the Gegenschein brightness) at a 



Galactic latitude of 50 ° (ILeiiiert et al.Ml998l). The inten sity of the airglow was subtracted us- 
ing Barbier's function ( Barbierl 1944 : Barbier et al. 1961 ). whereas the intensity of the multiple- 



scattering component in Earth's troposphere was not subtracted because of its weak contribu- 
tion in our data around the Gegenschein. The extinction coefficients of Earth's atmosphere were 
obtained through stellar photometry, while the effective optical depths for the extended sources 
were obt ained using the QDM_s ca code by adjusting the particle albedo u and the asymmetry 



factor Q (IHong et al 



1998, 



2OO2I) to match the brightness of previous zodiacal light observations 



(jLevasseur-Regourd &: Dumontlll98d ). Some of unknown parameters such as the intensity of the 
airglow at zenith and scaling factor of starlight brightness by Pio neer to our system were deter- 
mined empirically to fit our data to the zodiacal light brightness in iLevasseur-Regourd &: Dumont 
(jl98d ). We derived the Gegenschein brightness in the range of 165-190 Sqiq. Throughout this 
data reduction, we found that there is an uncertainty of ±10-15 S01O in the absolute brightness. 
This uncertainty is mainly caused by inadequate subtraction of atmospheric diffuse light sources 
and imperfect correction of the zenith extinction, although instrumental effects such as scattered 
light inside the camera cannot be ruled out. It is, however, emphasized that the relative intensity 
distribution of <20° structures is reliable in our data set because it cannot be influenced by the 
inevitable atmospheric and instrumental effects. 



3. Results 

3.1. Observational Results 

Figure O shows snapshots of the Gegenschein in ecliptic coordinates. The long axis was nearly 
parallel to the ecliptic and stretched in the east-west direction. The morphology was symmetrical 
with respect to the line at an ecliptic longitude with respect to the Sun of A — A0 = 180°, where A 
and A0 denote the ecliptic longitude and ecliptic longitude of the Sun. Looking at the morphologies 
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closely, we see a truncated oval shape whose eastern and western edges are angular rather than 
round. As the antisolar point is approached, the shape changes from truncated oval to round. 
There is a prominent peak at the antisolar point. 

It is clear that the brightness distribution shifted northward on 17 February and southward 
on 12 October (Figure [6j). We found that the Gegenschein morphology deviated to the south in 
August-October and to the north in February- April. This trend is consistent with the result of 
Ishiguro et al.l (|l998l ). who mentioned that the overall brightness distribution shifted southward in 
September but shifted northward in March. The annual variation can be explained by the three- 
dimensional zodiacal cloud model having a plane of symmetry deviated southward during August 



and N ovember and northward during February and April, as previously noticed in llshiguro et al, 
(| 19981 ). We discuss the annual variation in Section 3.2. 



Figure [7] shows the surface brightness profile of the Gegenschein on 12 October at three different 
longitudes. These profiles show a convex feature over the ecliptic latitude range /3 = —3° to 
1°. The flat-topp ed profiles were detected for the first time through the infrared observation of 



zodiacal emission (ILow et al.l 11984 ). They were also recorded for the Gegenschein (jlshiguro et al 
1999al : iMukai et al.l |2003| ) . Recent theoretical studies have showed that features in the profiles 



are ass ociated with materials gener ated by recent catastrophic collisions in the main-belt asteroid 



2003 



region (jNesvorny et al. 

A — A0 = 180°. Instead, it exhibits a steep increase at /3 



2009 ). However, the convex structure is not clear in the profile of 

0°. 



Despite the annual variation in the overall shape, we found that the maximum position coin- 
cided with the antisolar point through the observation runs. Neither east-west nor north-south 
displacement was detected in the profiles at A — A© = 180°. Figure [8] shows the ecliptic latitude of 
the maximum int ensity at X — Xp) = 180°. For comparison, we show photoelectric and photographic 



data co mpiled by RoosenI (Il970bl') as well as p revious CCD d a ta. O ur res ults are incon s istent with 



those in iRoosenI (|l970bl ). Ijames et al.l (119971 ) , llshiguro et al.l (119981 ) , and 



Mukai et al 



torn but 



they are consistent with the recent measurement from space by SMEI (jBuffington et al.ll2009l ). It 
should be emphasized that our measurement is the most reliable of all observational data so far 
because of the better spatial resolution. 



3.2. Comparison with 3-D Empirical Dust Distribution Models 

In this chapter, we analyze the annual variation in the morphology and the peak position by 
conducting three-dimensional zodiacal cloud model simulations. 

The surface brightness of the zodiacal light as seen by an observer in the ecliptic plane at an 
ecliptic latitude f3 and an ecliptic longitude X — Xq with respect to the Sun can be described as 
a double integral over the size distribution and along the line of sight. As an approximation, we 
assume the following: 
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1. The size distribution of dust particles in the zodiacal cloud is independent of their position 
in interplanetary space. 

2. The albedo and scattering phase functions of dust particles are independent of position. 

3. The zodiacal cloud has a plane of symmetry and is axisymmetric with respect to an axis 
perpendicular to the plane of symmetry. 

4. Dust grains are sufficiently far from the Sun that we can regard the Sun as a point source. 



Under these assumptions, the brightness of the zodiacal light is given by 

poo /'OO 

/zl(A - A0,/3) = / / F{rh)n{x,y,z) f{s)asca{s)<l){9;s)dsdl, 



JO 



(1) 



where dl and 6 denote the line element and the scattering angle, respectively, and s is the particle 
radius, /(s) is the differential size distribution of the particles with a radius of s. CTscais) is the 
cross section of the particles. 4>{6; s) characterizes the distribution of scattered intensity with 9. 
The incident solar flux F{rh) can be replaced by FQ{rQ/rh)'^, where Fq denotes the solar flux at 
the heliocentric distance rh = tq. n(x,y,z) is the number density of the cloud particles in ecliptic 
coordinates with respect to the Sun, where the x axis points from the Sun toward the vernal 
equinox, the z axis is perpendicular to the ecliptic plane, and the y axis completes a right-handed 
orthogonal coordinate system. The geometry and the notations are summarized in Figure [9l 

For simpliflcation, the mean volume scattering phase function ^{0) and the mean total scat- 
tering cross section o-gcat are defined by 



/(s) CTscatis) (t){9; s) ds = (Tscat ^(9), 



(2) 



with a normalization condition, that is. 



<l>{9)dn = 1. 



(3) 



47r 



The me a,n volume scatteririg phase function was empirically d erived from the zodiacal light obser- 
vations (jLeinert et al.lll976l : lHong]|l985l : iLamv &: PerrinI 119861). Figur e [TOl (left) shows an exampl e 
of the mean volume scattering phase function given by iHond (jl985l ) and iLamv &: PerrinI (jl986l ) . 
It has a strong peak in the forward direction, an isotropic component at intermediate scattering 
angle, and a slight enhancement in the backward direction. 

The spatial distribu tion of dust particles has been discussed (jGiese et al ] ll98d :l iGiese &: Kneissell 



19891 : iKelsall et al.lll998l ). In the early years of research, the spatial distribution was of a form sep- 
arable into radial and vertical terms, that is. 



no 



- ) 



(4) 
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where no is the reference dust number density on the symmetry plane at hehocentric distance 
''h = ^0- The radial power law is motivated by the radial distribution expected for particles under 
the effect of Poynting-Robertson drag, which results in = 1 for dust particles in circular orbit; 

describes how the dust number density reduces as we move away from the symmetry plane. 
Various axisymmetric distribution models were proposed. Among them, the ellipsoid model. 



h{l3') = [(1 + 6.5 sin /3' 



./\21 



-0.65 



and the fan model, 



h{f3') = exp(-2.1|sin/5'|), 



(5) 



(6) 



obtained igym ~ 



2° and 



sym 



80° , and iMukai et al.l (j2003l ) , who obtained 1 



were often applied in phenomenological model calculation in the literatu re (e.g.. iGiese &: Kneissel 
1983 ). The plane of symmetry was derive d by different authors (see, e.g.. iKwon et al.l (|2004l ). who 



"-sym 



2.03° ± 0.50° 



and Ogym = 54°-64°). A sophisticated model, in which localized dust condensations as well as the 
axisymmetric smooth component are considered, was suggested on the basis of infrared observations 
with the Cosmic Ba ckground Explorer's Diffuse Infrared Background Experiment (COBE/DIRBE) 
(jKelsall et al.lll998l ). It consists of a smooth diffuse component that has a plane of symmetry of 
^sym = 2.03° lb 0.017° and Ogym = 77.7° it 0.6°, three asteroidal dust components, and circumsolar 
ring components. 



We began with a comparison of the ellipsoid model and the fa n model. We assumed a plane of 



symm etry of i 



sym 



2° and^^, 



sym 



80° based on previous studies (jKwon et al.ll2004l : lKelsall et al 



19981 ). The radial power law exponent u = 1.3 is assumed. We confirmed that u does not change 
the Gegenschein morphology signifi c antly within ~30° from the antisolar point. The mean volume 
scattering phase function in iHone) (|l985l ) was applied. We should note here that the pow er-law 
expon ent i' is insensitive to the morphology of the Gegenschein, as already pointed out in iHong 
(|l985l ). For the reason, we did not give much attention to i' in this paper. The results are shown 
in Figure [TTl The systematic change in north-south asymmetry is explained qualitatively by these 
models. The model morphologies are, however, different from the observed ones in that they have 
protruding edges toward the eastern and western directions (see the arrows in Figure [TT]) . No 
peak in the antisolar direction appears in these modeled images. These mismatches cannot be 
explained by different values of Zgym and r^sym, suggesting that the classical models are inadequate 
for explaining the observed distribution of interplanetary dust particles. 



Next, we applied the infrared model of lKelsall et al.l (|l998l ). Because the model was constructed 
to fit the data in near- and far- infrared wavelengths, no scattering phase functi o n was provided for 
optical wavelengths. We again applied the scattering phase function in iHongi (|l985l ). The upper 
panels in Figure [T2] show the results using only the smooth component of the infrared model, while 
the middle panels show the results using all the components. From the comparison between top and 
middle panels, it is clear that the angular shape is associated with fine-scale structure, in particular, 
the dust-band components. Both the north-south asymmetry and the angular edge along the east- 
west direction were reproduced by the model (Figures 112( c) and (d)). It is interesting to note that 



-li- 



the model was constructed based on infrared observations and yet fit the optical zodiacal light as 
well. In addition, the infrared observations covered the zodiacal emission at solar elongation of 64°- 
124°. Both the observed wavelengths and the solar elongation angle of the infrared observations 
differ from those of our observations. Nevertheless, the infrared model matches the Gegenschein in 
the visible wavelength range. 

In Figures [12] and [6l we would like to draw attention to the antisolar point, where the brightness 
spike is clearly seen in our observational images but not visible in the model image. The difference 
is also seen in Figure \T3\ where the observed brightness profile around the antisolar point does not 
match the observed one. This is probably because the backscattering part of the scattering phase 
function is inadequate for reproducing the observed spikes. Since the scattering phase function 
has so far been derived using low-resolution data, it is likely that the backscattering part within 
~1° has not been considered y et. For this reason , we m odified the scattering phase function by 



multiplying it by a term used in iKaasalainen et al.l (|2003l ) , 



^\e) = {l + A exp(-^) ^ m, (7) 



where A and B are fitting parameters and ^(9) is the original scattering phase function in iHong 



(| 19851 ). We derived ^ = 0.15 and B = 0.021 (or 1.2°) to match the observed images. With the 
modified scattering phase function, we fit the surface brightness profiles (see Figure [12] (e)-(f) 
indicated by arrows, and also Figure [T3|) . 



4. Discussion 

The question we have to address here is whether the brightest spot on the Gegenschein was 
caused by scattered light from interplanetary dust particles, by an instrumental effect (i.e., uneven 
sensitivity of CCD pixels or vignetting), or by other light sources such as stars and galaxies. Since 
we detected the brightness enhancement in different positions on the CCD, it is unlikely that the 
brightness enhancement was caused by instrumental effects. In addition, because the antisolar 
position moves relative to stars and galaxies at a rate of l°/day, we can distinguish zodiacal light 
components from stellar objects by comparison with images taken on different nights. For these 
reasons, it is clear that the brightness enhancement originates from the interplanetary dust cloud, 
because it always appears at the antisolar point. 

Our observational result indicates that the volume scattering phase function sharply increases 
close to a phase angle of zero (that is, phase angle = 180°- scattering angle 9=0°). The sharp 
increase (hereafter referred to as the opposition spike) is a common phenomenon inherent to airless 
bodies in the Solar System, such as the Moon, asteroids, satellites of the major planets, Saturn's 
rings (Mishchenko 1993), and Kuiper Belt objects (Sheppard and Jewitt 2002); it is characterized by 
a nonlinear increase in surface brightness as the phase angle decreases to zero. Through this effect, 
the amplitude and the width of the scattered light change significantly within l°-5°. There are 
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two main physical mechanisms to explain the opposition effect on a surface covered with regolith: 
shadow-hiding and coherent backscattering enhancement. Shadows directly opposite the Sun are 
hidden by the scatterers, and the absence of visible shadow makes the area brighter on average 
(the shadow-hiding effect). On the other hand, coherent backscattering is an effect that describes 
the appearance of an intensity cone when photons traveling in a time-reversed path self-interfere 
constructively in the backscattered direction. 

The opposition effect has been detected to occur on airless bodies in the Solar System. Shadow- 
hiding and coherent backscatteri ng enhancement have been proposed as reasons for the opposition 



effect observed in regolith layers (jMuinonen et al.ll2002l ) . It is interesting to note that the opposition 



effect was detected in the interplanetary dust cloud from our observations. This evidence may sug- 
gest that scattering by a single interplanetary dust particle induces an opposition surge. Shkuratov 
& Helfenstein (2001) predicted from their model, in which the aggregate structure of regolith is con- 
sidered, that a single regolith particle could exhibit its own opposition effect. Hapke (200 2) also em- 



phasiz ed that multiple scatterings can occur between different parts of a single particle. iMin et al 



(|2010l ) showed that the large (> 100 /im) dust grains covered by regolith p articles are responsr 



bleto the strong backward scattering observed in the Fomalhaut debris disk (jLe Bouquin and al. 



2003 ). Further researches involving multiband observations, polarimetric observation and theoret- 
ical studies are required to study the mechanisms responsible for the opposition effect observed in 
the Gegenschein. Although we cannot conclude whether the opposition effect in the Gegenschein 
is caused by the shadow-hiding effect or interference or both, we can safely conclude that the effect 
is induced by scatterers whose sizes are significantly larger than the optical wavelength. 

The albedo of interp laneta r y dus t particles leads us to consider their origin. This has been 



investigated in literature; ICookl (|l978l ) studied the albedo using micrometeoroid experiment data 
onboard Helios A and Pioneer 10 as well as zodiacal light observational data and derived an albedo 
ranging from 0.07 at 0.1 AU from the Sun through 0.006 at 1 AU and down to 0.0001 at 3.3 
AU. Later, Hanner (1980) reanalyzed the data to obtain a mean geometric albedo of 0.24 for dust 
particles near 1 AU. Dumont & Levasseur-Regourd (1988) and Renard et al. (1995) studied zodiacal 
light and zodiacal emission by using the node of lesser uncertainty method, and they obtained an 
albedo at 1 AU of 0.08-0.15 at a scattering phase function of 90°, which translates into a geometric 
albedo of 0.16-0.30 with the scattering phase function by Hong (1985). Lumme & Bowell (1985) 
derived an albedo of 0.04 using an empirical polarimetric method. All these measurements are 
inconsistent with each other and have not yet led to a discussion about the origin. 

Our observational data can be used to derive an accurate albedo value in combination with 
the infrared model. Because our observed quantity is obtained at the opposition direction from the 
Sun, it gives straightforward input for deriving the geometric albedo of interplanetary dust particles 
around the Earth orbit (Appendix A.) because it is defined at zero phase angle. In addition, the 
good consistency between the optical and infrared spatial distributions of dust particles suggests 
that the effective size of dust grains in the optical wavelength region is almost identical to that at 
infrared wavelengths. This result supports the idea that large grains (>10 ^m) are responsible for 
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the scattering cross section of zodiacal light (jGiese &: Gruenlll976l ) 



In light of this large-particle hypothesis, let us derive the geometric albedo. We rewrite Eq. 



(1) at the antisolar point using the geometric albedo ^j,(180°) (jHanner et al.lll98lh as follows: 

/zl(180°,0°)= r F{rMx,y,z)^^^^^dl. (8) 
Jo ^ 

Because n(x, y, z) is given by the DIRBE infrared model while the intensity of zodiacal light at 
opposition is obtained by our observation (175±15 Sqiq), we can deduce a mean geometric albedo of 
^p(180°) = 0.07 ± 0.01. The DIRBE infrared model also suggests that about 7% of the geometric 
cross-sectional area accounts for asteroidal particles, particularly the inner dust ba n d from the 



Kariri collisional family (whose mean geometric albedo is 0.20-0.22; [Harris et al.l (j2009l ) : lUsui et al 



(j201ll )) and the circumsolar resonance ring possibly originating from the asteroidal belts. Given 



that 7% of the geometric cross section is accounted for by particles with ^p(180°) = 0.20-0.22, we 
found that the remaining zodiacal light (the smooth component in the DIRBE model) should have 
a geometric albedo of 0.06 ± 0.01. 

Our result is consistent with the fact that the micrometeorites accumulated around Earth or- 



bit have low albedo (<0.1: iBrownled (jl978l ): iHannerl (jl980l )). Among Solar System minor objects, 
comet nuclei and C- and D-type asteroids show low a lbedo, that is, .02-0.06 for comet nuclei 
(jLamv et al.ll2004l ) and 0.07 ±0.04 for C-type asteroids (jUsui et al.ll201ll ). Recent dynamical simu- 
lations of interplanetary dust particles suggest that 85%-95% of the observed mid-infrared emission 
i s produced by partic le s from Jupiter-fain ily comets and <10% by dust from long-period comets 
(jNesvorny et al.ll2010l ). Ijewitt et all (j2013l ) conclude that <3% of the interplaneta ry dust particles 



are produced by events l ike P/2010 A2 pro bably caused by breaku p of an asteroid (jSnodgrass et al 



2010 



Jewitt et al. 



2011 



On the contrary, iTsumura et al.l (120101 ) insist that high-albedo particles 
dominate near Earth orbit based on their spectroscopic observation at near- infrared wavelengths. 
Although further investigation is required on the fraction of asteroidal particles with respect to 
cometary particles, it seems plausible that majority of interplanetary dust particles have a compo- 
sition similar to carbonaceous chondrites rather than siliceous particles. 



5. Summary 

In this paper, we outlined the instrumentation, the data reduction, and the observations of the 
Gegenschein using the hand-made observation system WIZARD. Our main findings are as follows: 

• A narrow brightness enhancement was superimposed on the smooth gradient of the Gegen- 
schein at the antisolar point. Hence, the position of maximum brightness exists at the antisolar 
point throughout the year. 

• The Gegenschein morphology changed according to the orbital motion of the Earth. It is well 
produced by a model constructed using infrared observations by COBE/DIRBE. 
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We hypothesized that the zero phase angle spike is caused by coherent backscattering and/or 
shadow-hiding effects on the rough surfaces of individual dust particles. The above results indicate 
that the effective sizes of interplanetary dust particles responsible for optical zodiacal light and 
infrared zodiacal emission are much larger than those corresponding to the observed wavelengths. 

Finally, we derived the geometric albedo of interplanetary dust particles using our observational 
result and found that it is 0.06 it 0.01 for the smooth component of the zodiacal cloud. The derived 
albedo is in accordance with that of the dark micrometeorites accumulated around earth orbit and 
the observed cometary dust particles. We concluded that chondritic particles are dominant near 
Earth space. 
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A. Appendix A: Line of Sight Depth of Zodiacal Light and Zodiacal Emission 

Because the intensities of optical zodiacal light and infrared zodiacal emission are quantities 
integrated along the line of sight, it is difficult to specify the particles we observe. We consider the 
contribution of the elements at different distance from the observers to the total flux of optical zodi- 
acal light and infrared zodiacal emission. Figure [HI show the contribution in different wavelengths 
at A — A0 = 90° and A — A© = 180° on the ecliptic plane. We found that the line of sight depth 
becomes shallowest at 12/im because hot dust particles (~280 K) near Earth orbit are efficient 
infrared emitters. On the other hand, distant particles contribute well at longer wavelength. The 
depth in the optical wavelength is similar to that at 25/.fm. It seems that we derived the albedo of 
interplanetary dust particle within ~ 2 — 3 AU from the Earth. 
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Fig. 1. — Photograph {top) and schematic view (bottom) of WIZARD. It consists of a cryogenic 
dewar and optics. These two units are connected with a behows tube and mounted on a thick 
honeycomb plate. The mechanical shutter and glass filter are placed at the pupil position. The 
system is mounted on an equatorial mounting for star-tracking. 
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Fig. 2. — Transmittance of the filter we applied for the Gegenschein and zodiacal light observations. 
This filter is designed to avoid the prominent natural emissions and artificial emissions, indicated 
by solid lines and dashed lines, respectively. They consist of O I (557.7, 630.0, and 636.4 nm), Nal 
(589.0 and 589.6 nm), and Ha (656.3 nm) for the natural emissions and Hgl (435.8, 546.1, and 
579.1 nm) and Nal (589.0 and 589.6 nm) for the artificial emissions, respectively. 
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Fig. 3. — Vignetting factor of WIZARD no rmali zed at the optical c enter (thick solid line). Two vi- 
gnetting functions in llshiguro et al.l ()1999al ) and I James et al.l (j 19971 ) are also drawn for comparison. 
It is clear that WIZARD optics has little vignetting within 20° from the center. 
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Fig. 4. — Ecliptic longitude (top) and Galactic latitude (bottom) of the antisolar point as a function 
of the day of the year. Crosses denote the ecliptic longitude and Galactic latitude when we made 
the observations. The antisolar point lies on the Galactic plane on June 21 and December 22. 
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Fig. 5. — Comparison between the instrumental magnitude defined as — 2.51og(DN) + constant 
and color-corrected catalog magnitude of field stars, where DN denotes CCD count. Since these 
stars exist near the zenith (within an airmass of 1.07), we did not correct for the zenith extinction 
term (see Section [2.2p . 
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Fig. 7. — Surface brightness profile along the ecliptic latitude at A — A© = 180° (top), 170° (middle), 
and 160° (bottom). For clarity, the middle and lower profiles are shifted by 10 S01O and 20 Sqiq, 
respectively. 
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Fig. 8. — Observed latitude of the maximum intensity of the Gegenschein. In the top panel, the 



( Roosen 


1970b: 


Ishieuro et al. 


1998) 



compa riso n, we show the pos ition predicted using a plausible model on the basis of iMukai et al 
(j2003l ) and lKwon et al.l (j2004l ) (solid line). In the model, t he pla r ie of s ymmetry for isym = 2° and 



''sym 

assuming i/ 



80° is considered. The scattering phase function in iHongl (119851 ) was applied for the model, 
1.3. 
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Fig. 9. — The viewing geometry for light scattered by a dust volume element. The scattering angle 
9 is measured from the line of sight to the antisolar direction at the volume element. 
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Fig. 10. — (Left) The scattering phase functions from lHona (j 19851 ) and lLamv &: PerrinI (j 19861 ) for a 
radial power law v = \. [Right) Comparison between scattering phase functions obtained by other 
authors and us. 
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Fig. 11. — An isophoto map, in arbitrary units, of the Gegenschein at two different days of year 
based on th e elhpsoid model (top) and the fan model (bottom). The scattering phase function of 
Hong) (jl985l ) was applied. 



- 29 - 



Kelsall model smooth component 
(without opposition surge) 



KeSsal! mode! smoolh component 
(without opposition surge) 




Kelsall mode! (without opposition surge) Kelsall model (without apposition surge) 




February 1 7 j 



■ ■ ■ — ■ ■ ■ 



■ I — -I. 



l—-^-:r:.r-T.Tr'" October 12 



Kelsall model (with apposition surge) 



Kelsall model (with opposition surge) 



200 190 180 170 160 1502 

X-k.. 



October 1 2 



20 
10 

P 

'10 

-20 
20 




200 190 180 170 160 150 



Fig. 12. — An isophoto map, in arbitr ary units, of the G egenschein at days from two different sea- 
sons based on the infrared 3-D model (iKelsall et al.lll998l ). We apphed only the smooth component 
in the top panels and all component s including d ust bands in the middle and bottom panels. The 
original scattering phase function of iHongl (|l985l ) was applied in the top and middle panels, while 
the modified scattering phase function was applied in the bottom two panels. 
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Fig. 13. — Comparison between the observed brightness profiles w ith those obt ained from the model 
(solid lines). We applied the original scattering phase func tion of iHongj (|l985l ) assuming v = 1.3 in 
the left panel and a modified scattering phase function of iHongj (j 19851 ) with 15% enhancement at 
opposition in the right panel. 
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Fig. 14. — The line of sight depth of optical zodiacal light and infrared zodiacal emission at at 
A — A0 = 90° and A — A© = 180° on the ecliptic plane. 
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Table 1: Summary of the WIZARD performance 



Field-of-view 


49°x98° 


Pixel resolution 


1.435' 


Gain factor 


0.825 (e~/ADU) 


Readout noise 


20-25 (e") 


Stability of the CCD zero level 


<10 (e-) or 

<0.6% of the Gegenschein brightness 


Vignetting 


30 % (at 47° from the optical center) 


Pixel-to-pixel responsivity variation 


<0.3% 
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Table 2. Observation summary 



Date 


UT" 






c 


z 


d 




bf 


March 01 


10: 


:05 





.991 


160, 


,36 


14. 


,04 


241 


+55 




10: 


;29 





.991 


160, 


,37 


12, 


,21 


241 


+55 


March 02 


09: 


:35 





.991 


161, 


,34 


18. 


,99 


242 


+55 




09: 


:58 





.991 


161, 


,36 


15. 


,25 


242 


+55 




10: 


:45 





.991 


161, 


,39 


12. 


,82 


242 


+55 




11: 


:35 





.991 


161, 


,42 


19. 


,41 


242 


+55 


August 25 


09: 


:58 


1 


.011 


331, 


.75 


31. 


,33 


49 


-50 




10: 


:26 


1 


.011 


331, 


,77 


30. 


,68 


49 


-50 




10: 


:49 


1 


.011 


331, 


,78 


31. 


,28 


49 


-50 


August 27 


10: 


:06 


1 


.010 


333, 


,69 


30. 


,28 


52 


-51 




10: 


:29 


1 


.010 


333, 


,70 


30. 


,01 


52 


-51 


February 15 


09: 


:41 





.988 


145, 


,98 


14. 


,87 


223 


+46 




10: 


:00 





.988 


145, 


.99 


11. 


,06 


223 


+46 


February 16 


08: 


:23 





.988 


146, 


,94 


32. 


,64 


224 


+47 




09: 


:01 





.988 


146, 


,96 


23. 


,86 


224 


+47 




09: 


:41 





.988 


146, 


,99 


15. 


,03 


224 


+47 




10: 


:05 





.988 


147, 


,01 


10. 


,39 


224 


+47 


February 17 


09: 


:09 





.988 


147, 


,98 


22. 


,17 


225 


+47 




09: 


:27 





.988 


147, 


,99 


18. 


,18 


225 


+47 




09: 


:44 





.988 


148, 


,00 


14. 


,58 


225 


+47 


April 21 


09: 


:28 


1 


.005 


211, 


,55 


34. 


,34 


327 


+48 




10: 


:05 


1 


.005 


211, 


,58 


32. 


,07 


327 


+48 




10: 


:24 


1 


.005 


211, 


,59 


31. 


,87 


327 


+48 




10: 


:42 


1 


.005 


211, 


,60 


32. 


,30 


327 


+48 


April 22 


09: 


:24 


1 


.005 


212, 


,52 


35. 


,01 


328 


+47 




09: 


:42 


1 


.005 


212, 


,54 


33. 


,52 


328 


+47 


September 12 


11: 


:53 


1 


.006 


349, 


,98 


33. 


,37 


77 


-59 




12: 


:12 


1 


.006 


350, 


,00 


36. 


,79 


77 


-59 


September 13 


11: 


:25 


1 


.006 


350, 


,94 


28. 


,71 


79 


-59 




11: 


:45 


1 


.006 


350, 


.95 


31. 


,80 


79 


-59 




12: 


:04 


1 


.006 


350, 


.96 


35. 


,12 


79 


-59 


October 13 


08: 


:52 





.998 


20, 


.29 


21. 


,87 


133 


-54 




09: 


:13 





.998 


20, 


.31 


17. 


,82 


133 


-54 
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Table 2 — Continued 



Date 
















October 14 


09:22 


n nn'7 

0.997 


o 1 on 

21.30 


15.95 


134 


—54 






09:41 


n nn*? 

0.997 


m on 

21.62 


1 o 1 n 

13.19 


134 


—54 






09:57 


n nn*? 

0.997 


n 1 o o 

21.33 


11.79 


134 


—54 






10:13 


0.997 


21.34 


11.58 


135 


-54 


November 


11 


09:59 


0.990 


49.29 


2.77 


163 


-34 






10:18 


0.990 


49.30 


3.70 


163 


-34 






10:36 


0.990 


49.31 


7.53 


163 


-34 


November 


12 


10:16 


0.990 


50.30 


3.14 


164 


-34 






10:33 


0.990 


50.32 


6.74 


164 


-34 






10:49 


0.990 


50.33 


10.41 


164 


-34 


November 


13 


10:19 


0.989 


51.31 


3.54 


165 


-33 






10:52 


0.989 


51.34 


11.02 


165 


-33 


2006 November 


18 


08:47 


0.988 


55.77 


18.76 


168 


-29 






09:05 


0.988 


55.78 


14.53 


168 


-29 


November 


23 


08:41 


0.987 


60.81 


20.38 


171 


-25 






08:59 


0.987 


60.83 


16.16 


171 


-25 


November 


27 


12:01 


0.987 


65.00 


26.21 


173 


-22 



Note. — "Median time of exposure. 
''Heliocentric distance of the Earth in AU. 

'^Ecliptic longitude of the antisolar point with the equinox of J2000 
in degrees. 

'^Zenith distance of the antisolar point in degrees. 
'^Galactic longitude of the antisolar point in degrees. 
■'^ Galactic latitude of the antisolar point in degrees. 



